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ABSTRACT 



Context. Relating in situ measurements of relativistic solar particles to their parent activity in the corona requires understanding the 
magnetic structures that guide them from their acceleration site to the Earth. Relativistic particle events are observed at times of high 
solar activity, when transient magnetic structures such as Interplanetary Coronal Mass Ejections (ICMEs) often shape the interplan- 
etary magnetic field (IMF). They may introduce interplanetary paths that are longer than nominal, and magnetic connections rooted 
far from the nominal Parker spiral. 

Aims. We present a detailed study of the IMF configurations during ten relativistic solar particle events of the 23rd activity cycle to 
elucidate the actual IMF configuration guiding the particles to Earth, where they are measured by neutron monitors. 
Methods. We use magnetic field (MAG) and plasma parameter measurements (SWEPAM) from ACE, and determine interplanetary 
path lengths of energetic particles through a modified version of the velocity dispersion analysis based on energetic particle measure- 
ments with SoHO/ERNE. 

Results. We find that the majority (7/10) of the events is detected in the vicinity of an ICME. Their interplanetary path lengths are 
found to be longer (1.5-2.6 AU) than those of the two events propagating in the slow solar wind (1.3 AU). The largest apparent path 
length is found in an event within the fast solar wind, probably due to enhanced pitch angle scattering. The derived path lengths imply 
that the first energetic and relativistic protons are released at the Sun at the same time as electron beams emitting type III radio bursts. 
Conclusions. The timing of the first high-energy particle arrival at Earth is dominantly determined by the type of IMF in which the 
particles propagate. Initial arrival times are as expected from Parker's model in the slow solar wind, and significantly larger in or near 
transient structures such as ICMEs. 

Key words. The Sun: solar-terrestrial relations - The Sun: heliosphere - Method : data analysis 



^ 1. Introduction 



00 

d 



> 
X 



The solar activity has consequences for the entire inner helio- 
sphere. Besides high-energy photons, one can distinguish two 
energetic phenomena directly affecting the teiTestrial environ- 
ment: coronal mass ejections (CMEs), disturbing the magne- 
tized environment, and the energetic particle events, impacting 
the Earth's atmosphere and affecting the ionized environment. 
The most energetic particles that the Sun produces are relativis- 
tic protons up to about 10 GeV. At these relativistic energies, so- 
lar particles penetrate the magnetosphere and impact the Earth's 
atmosphere. The atmospheric interaction of the GeV particles 
produces secondary particles through nuclear cascades that are 
detected at the ground level by neutron monitors (NMs). Hence, 
the name of ground level enhancement or GLE. Only 70 events 
have been repotted since 1942, and have been largely studied in 
order to address the acceleration and propagation of these parti- 
cles from Sun to the Earth. 

Solar energetic phenomena may accelerate particles up to 
relativistic energies through the coronal shock driven by the 
CME (Vainio & Laitinen| [20071 |Sandroos & Vainio| [2009ll, or 
through magnetic reconnection during a flare ([Arzner & Vlahos 
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IDmitruk et al] [20041 [Turkmani et al.| |20051 [Drake "eTaL 



Energetic particles are guided by the interplanetary magnetic 
flux tube connecting the acceleration site to the Earth. Most pre- 
vious studies are based on the assumption that energetic par- 
ticles propagate along the nominal Parker spiral. Nevertheless, 
two main problems arise under this assumption. First, it implies 
that the parent active regions are located near 30° - 80° West. 
However, observations show that the distribution in longitude of 
the parent active regions is broad, ranging from 90° East to more 
than 120° West ( |Cliver et al.[[T9Ml[Kahler etaD T984). Second, 
the time measured between the radiative signatures of particle 
acceleration in the Sun's corona and the in-situ measurements 
is longer than it should be for a propagation in the Parker spiral 
(Debrun neretal.| ri997; Kahler et al., 2003; Tylka et al. 2003). 

Several scenarios have been proposed to explain both the 
delay and the connection problem during relativistic particle 
events. Indeed, the delay may be due to a late acceleration and/or 



injection phase during the flare ( Debrunner et al.[ 1997 Klein 
l et al.| [1999)1 or to a time extended acceleration at the bow shock 
of the CME (Reames, 1999; Gopalswamy, 20051. Since shock 
acceleration should inject particles over a large angular width, 
explaining also the connection problem, several studies favored 
the CME shock acceleration rather than the flare acceleration 
dCliverl [T9821 [Hudson et al.l|1982, Cliver et al. 2004). 

However, it has been known since early investigations of rel- 
ativistic solar energetic particle events that they often occur dur- 
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ing depressions of the galactic cosmic ray intensity. Such a de- 
pression is typically associated with a transient interplanetary 
magnetic field structure (see Carmichael 1962 and references 
therein), now known as an interplanetary coronal mass ejection 
(ICME), which is the interplanetary counter part of the CME 



observed in the solar corona ( Wimmer-Schweingruber et al.| 
|2006 '). While observed at 1 AU, ICMEs can still be connected 
at the Sun by one or two footpoints (e.g. |Crooker et al. |2008] l. 
Thus, relativistic particles could be detected at Earth in ICMEs 
when the parent active regions is far from the nominally well- 
connected western solar hemisphere ( jPebrunner et al. 1988| l. 
Detailed studies of individual events showed indeed that ener- 



getic particles can propagate within ICMEs (Richardson et al. 
[T99Tl|Larson"etaLl|T997i|Torsti et al.||2004| . 



Moreover, the specific magnetic topology of ICMEs mod- 
ifies the interplanetary propagation of energetic particles such 
as the directivity of the particles flux ( Krittinatham & Ruffolo 
2009) 1. This could explain peculiarities of directional distribu 



tions of relativistic protons observed at Earth ( Bieber et al. 
2002"; Miros hnichenko et"aLi [2005i [Ruff^olo et al.[ "'20%7'Sliz 
et al., 2008 ). Also, long geometrical paths in the interplanetary 



space inferred from velocity dispersion analysis ( Larson et al. 
1 1997 ') or from detailed timing comparisons ( [Masson et al.[|200^ 
point to particle propagation in non-nominal interplanetary mag- 
netic fields. 

Thus, the relationship between solar energetic particles mea- 
sured near 1 AU and the parent activity in the corona depends on 
a combination of processes: particle acceleration in the corona, 
access of the accelerated particles to open field lines, and propa- 
gation through interplanetary space along various possible mag- 
netic field configurations. 

In this paper we explore the interplanetary magnetic struc- 
tures and their impact on particle propagation during the most 
energetic solar particle events, the GLEs. Our analysis is based 
on two independent methods: the identification of the interplan- 
etary magnetic structures using magnetic field and plasma pa- 
rameters measured aboard the ACE spacecraft (Section |2]) and 
the velocity dispersion analysis of the initial proton arrival times 
based on ERNE/SoHO and neutron monitor data (Section [3]l. 
This is completed by comparing the deduced solar release time 
of protons with the time interval of electrons as deduced from 
type III bursts observations. Then we summarize and conclude 
in Section m 



2. The interplanetary magnetic structures guiding 
relativistic particles 

2. 1 . Characteristics of interplanetary structures 

The magnetic field vector at the spacecraft is conveniently char- 
acterized by the magnitude of the field, B, its orientation de- 
fined by the latitude, 9b, above the ecliptic plane and the longi- 
tude, (f>B, with respect to the Sun-Earth axis. In the case of the 
undisturbed solar wind, where the field line is a Parker Spiral, 
the magnetic field is weak (B ~ 5 nT at 1 AU) and strongly 
non coherent, with angles fluctuating around typical mean val- 
ues ^ 0° and (p^ - -45° or 135 ° (for an inward or outward 
sector, respectively, with GSE coordinates). When an ICME, ini- 
tially ejected from the solar corona, reaches the spacecraft, B 



increases strongly and it becomes more coherent (e.g., Burlaga 
1995[ l. The magnetic field orientation evolves also more coher- 
ently. 



In the solar wind, the temperature of protons is empirically 
dependent on the proton velocity (Lopez & Freeman} 1986 1. 
Elliott etal.| ( |2005| ) established the relation 



rexp[K] = 640 X Vsw[km-s"'] - 1.56 x lO^' 



(1) 



between the expected proton temperature, Texp, and the solar 
wind speed, Vsw- A deviation from this relationship allows us 
to identify a transient interplanetary structure. An ICME has 



typically a proton temperature lower than rexp/2 ( [Burlaga et al. 
TMTIIKlein & Burlaga] 1982} |Burlaga[T99T]), and this criterion is 



frequently used to define the extent of an ICME ( [Richardson & 
Cane 2010 and references therein). Moreover, if the magnetic 
field is also significantly more intense than in the solar wind, 
with a coherent rotation, this region defines a magnetic cloud. 

Differences between solar wind and ICMEs are also ex- 
pected in the proton /3 (ratio of proton over magnetic pressure, 
noted /3p). Typically /3p > 0.4 in the solar wind, while it is gen- 
erally smaller, typically in the range 0.01 < /3p < 0.4, within an 
ICME or a magnetic cloud, which are cooler and have stronger 
magnetic field. 

When a magnetic cloud (or an ICME) moves faster than the 
surrounding solar wind, a sheath is formed in front with en- 
hanced plasma density and magnetic field strength. Magnetic 
reconnection is typically expected between two magnetic struc- 
tures when different magnetic fields are pushed against each 
other. In situ evidences of such reconnection have recently been 



found in a magnetic cloud sheath (Chian & Mufioz 201 1 1. More 



generally, reconnection between the sheath and flux rope mag- 
netic field implies that the flux rope front can be progressively 
pealed, while at the flux rope rear the corresponding magnetic 
flux region remains (called a back region). Because of the differ- 
ent magnetic topologies, the back region is typically separated 
from the remnant flux rope by a current sheet. This back re- 
gion has intermediate properties between those of the solar wind 
and the magnetic cloud, since it originally belongs to the flux 
rope, but after reconnection it is also connected to the solar wind 
( [Dasso et al.[ [20061 [2007 ). 

A priori, solar energetic particles can travel in any of the 
interplanetary structures summarized above. So, depending on 
their injection site on the Sun, or their access to the structure, 
energetic particles could be observed in a parkerian solar wind 
or in an environment associated with an ICME (e.g., in the ICME 
sheath, inside a magnetic cloud flux rope, in the magnetic cloud 
back, or in the distorted solar wind that is perturbed by the fast 
transient). 

2.2. GLE selection and characterization of the interplanetary 
magnetic structures 

We initially select the relativistic events reported between 2000 
and 2006. We do not study the GLEs with a marginal increase 
of the neutron monitor count rate (< 5%). The list is thus re- 
strained to 10 events (Table [T]). The arrival times of the first de- 
tected protons at Earth have been already published by [Moraal[ 
& McCracken ( 201 1[ ), and we report their results in Table [IJ 

Om analysis is based on measurements of the plasma and 
magnetic field near the LI Lagrange point of the Sun-Earth 
system by the Advanced Composition Explorer (ACE) mission 
( Stone et al.[ 1998[): the magn etic data are obtained from the 
MAG instrument (Smit h et al.[[1998|l and th e plasma data from 
the SWEPAM instrument ( [McComas et al.| [1^98) . In order to 
compare with neutron monitor measurements at 1 AU, we have 
to propagate the ACE measurements to the Earth. We perform 
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Column 1 


2 


3 


4 


5 


6 


7 


8 


# GLE / date 


^onset(UT) 


Max. % 


Active region 


PS 


Vsw 


Connection Interplanetary 








increase 


location 


footpoints 


(km.s-') 




structure 


59/ 14 Jul. 


00 


10:31 


59 


N22,W07 


W43 


500 


poorly 


back region 


60/ 15 Apr. 


01* 


13:55 


237 


S20,W85 


W33 


640 


poorly 


disturbed solar wind 


61/ 18 Apr. 


01* 


02:36 


26 


behind limb 


W60 


360 


poorly 


ICME sheath 


63/ 26 Dec. 


01* 


05:40 


13 


N08, W54 


W51 


420 


well 


slow solar wind 


64 / 24 Aug . 


02* 


01:24 


14 


S02, W81 


W54 


400 


well 


slow solar wind 


65/28 Oct. 


03* 


11:12 


47 


SI 6, EOS 


W48 


450 


poorly 


back region 


66/29 Oct. 


03 


21:01 


35 


S15, W02 


W48 


450 


poorly 


magnetic cloud 


67/02 Nov. 


03* 


17:30 


39 


S14, W56 


W36 


600 


well 


back region 


69/20 Jan. 


05 


06:49 


5400 


N14, W55 


W39 


550 


well 


disturbed solar wind 


70/ 13 Dec. 


06* 


02:50 


92 


S06, W23 


W33 


660 


well 


fast solar wind 



Note : Results obtained from independent studies are synthesized in this table. Col. 1: the number and the date of the GLEs. Col. 2: the onset 
time of the first responding Neutron Monitor. Col. 3: Maximum percentage increase of the GLE ( [Belov et al.[|2010l >. Col. 4: the location of 
the parent active region. Col. 5: the solar longitude of the Parker spiral footpoint. Col. 6: The mean velocity of the solar wind computed on 
a half-day time interval for a quiet region of the IMF immediately preceding any transient magnetic structures transiting through the Earth 
during the GLE. Col. 7: The nature of connection between the parent active region and the Parker spiral footpoint. Col. 8: results of the 
identification of the magnetic topology of the interplanetary medium. The stars on the right of Col. 1 indicate the GLEs where interplanetary 
length and solar release time could be estimated (Table [2j. 



this time correction by assuming bodily propagation of the mag- 
netic field structures from LI to Earth at the average solar wind 
proton speed Vp measured with ACE / SWEPAM during the two 
hours preceding the arrival of protons at Earth. This implies a 
time shift At = Du/Vp where Du is the Ll-Earth distance. 

The magnetic field and the plasma parameters are drawn on 
adjacent panels, in order to compare their temporal evolution 
(e.g. Figs. 1 3 I. By over-plotting as a grey bar the 1 -hour-time 



interval starting with the first arrival time of relativistic particles 
at the first responding neutron monitor, we characterize the mag- 
netic structure in which the first relativistic protons arrive at the 
Earth. 

We localize their parent active region (col. 4 of Table [TJ and 
the longitude of the footpoint of the Earth-connected Parker spi- 
ral field line, given in col. 5 of Table [T] The longitude of the 
Parker spiral root on the solar source surface is at ( [Parker^ ^1961^ 



= 0(1 AU) + —(1 AU - Rs) , 



Vs 



(2) 



w 



where Q. - 2.6 x 10"^ rad.s ' is the angular speed of the Sun, 
Vsw the solar wind speed, assumed constant from the Sun to the 
Earth, /?, - 2.5Rq is the radius of the spherical source surface, 
and 0(1 AU) is the longitude of the spiral at 1 AU. We select 
<1)(1 AU) = as the reference longitude and O^ is positive west- 
ward from the Sun-Earth axis. Because of magnetic fluctuations, 
|Ippolito et al.| ( |2003] ) showed by a numerical study that the longi- 
tude of the Parker spiral footpoints at the source surface can not 
be estimated with an accuracy better than 6° - 10°. Moreover, 
coronal magnetic field extrapolations showed that an open coro- 
nal flux tube rooted in an active region can spread by several 
tens of degrees in longitude westward and eastward ( |KIein et alT 



|2008| l. Therefore, we consider that the parent active region is 
well-connected to the Earth by the Parker spiral if its longitude 
is within the range [<J>s - 30°; + 30°], otherwise we define 
the parent active region as being poorly connected (col. 7 of 
Table [TJ. The path lengths along the nominal Parker spirals in 
this sample range from 1.06 AU (Vsw - 660kms"') to 1.15 AU 
(400 km.s"'), i.e. they are essentially 1.1 AU for each event. 



2.3. Magnetic configurations guiding relativistic particles 

Our analysis leads us to separate the magnetic field configura- 
tions along which relativistic particles reached the Earth into 
three distinct subsets (col. 8 of Tab le [T): t he nominal solar wind 
with a Parker spiral field (Section 2.3. l|l, a transient magnetic 
structure related to an ICME (Section |2.3.2| i, and a third cate- 
gory, typically a highly disturbed solar wind in the vicinity of an 
ICME (Section [031 l. 



2.3.1. The Parker spiral 

In our sample, only 3 GLEs reach the Earth along a typical 
Parker spiral: Dec. 26, 2001 (GLE 63), Aug. 24, 2002 (GLE 64) 
and Dec. 13,2006 (GLE 70). 

Between 02:00 - 24:00 UT on Dec. 26, 2001, the magnitude 
of the magnetic field is fluctuating around 5 - 6 nT (Figure [TJ. 
The magnetic field orientation is typical of the Parker spiral 
(6b ~ 0° and 0b ~ -45°). Moreover, the expected temperature is 
almost the same as the observed one and 0.4 < /3p < 1. A small 
magnetic cloud could be present between 21:00 UT on Dec. 25 



and 02:00 UT on Dec. 26, but it was not listed by Richardson & 
Cane]( [20T0) l and it ends about 4 hours before the arrival of rel- 
ativistic particles. Thus, the relativistic particles impacting the 
Earth at 05:40 UT on Dec. 26, 2001 (GLE 63), propagate along 
the Earth-connected Parker field line. 

The two other events consistent with a propagation along 
the nominal Parker spiral are GLE 64 and GLE 70. The mag- 
netic field orientation, the plasma beta and the similarity of the 
expected and the actually observed temperature all clearly ar- 



gue in favor of this interpretation (see Sections A. 1 . 1 
Figures [ATT||A.2|i. 



A. 1.2 and 



For these 3 events, no ICMEs have been identified few days 
before the particle arrival (Richardson & Cane 2010[ l, confirm- 
ing that the interplanetary medium is not disturbed. Moreover, 
the proximity in longitude of the parent active region and the 
footpoint of the nominal Parker spiral is consistent with this re- 
suh (Table [1). 
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Fig. 1. Interplanetary data around GLE 63 (Dec. 26, 2001): ex- 
ample of a GLE present in a typical Parker spiral. From the top 
to the bottom, panels display respectively: the magnitude of the 
magnetic field (nT); its latitude 9^ and its longitude (p^ (given in 
GSE coordinate system); then on the same panel are over plot- 
ted the temperature of protons measured by the spacecraft {black 
dots) and the expected temperature {gray dots), and finally the fi 
of protons. An horizontal line is added for 9b - Q, (pB - -45° 
and jSproton - 0.4 to mark characteristic values. The grey-scale 
rectangle on the 5 panels corresponds to the 1-hour time interval 
starting at the first arrival time of relativistic protons at Earth. 
The time axis is as measured at Earth. ACE measurements at the 
LI Lagrangian point have been shifted by the appropriate travel 
time (see Section|Z2). 



2.3.2. ICME, magnetic cloud or back region 

For the five GLEs occurring on Jul. 14, 2000 (GLE 59), on Apr 
18, 2001 (GLE 61), on Oct. 28, 2003 (GLE 65), on Oct. 29, 
2003 (GLE 66), and on Nov. 2, 2003 (GLE 67), we find that 
particles propagate in a transient magnetic structure, such as in 
the sheath of an ICME, in a magnetic cloud or in the back region 
of a magnetic cloud. 

The GLE on Nov. 2, 2003 is detected at Earth at 17:30 UT 
(Figure |2]i. On that day, the magnetic field strength {B ^ 5 nT) 
and its longitude (fluctuating around 0^ - -50 °) are typical 
of a solar wind structured by the Parker spiral. However, 9b is 




0.01 L 
31-Oct 



01 -Nov 02-Nov 03-Nov 04-Nov 
Start Time : 03/10/31,0000:00 UT 



Fig. 2. Interplanetary data around GLE 67: it is located in the 
back region of a magnetic cloud. Same drawing convention as 
Figure [T] 



not aligned with the ecliptic plan (^ -50 °) as one expects for a 
Parker spiral. Moreover, both Tobs < T^exp and fip < 0.4 are in- 
consistent with a quiet solar wind. Indeed, an ICME is detected, 
starting on Oct. 31, 2003 at 2:00 UT and ending on Nov. 2, 2003 
at 00:00 UT 

Within the ICME, from ^ 8:00 UT on Oct. 31 to 1:00 UT on 
Nov. 1, both 9b and (pB have a relatively coherent global evolu- 
tion. Moreover, Tobs < Tsj^p/l, implying that a magnetic cloud is 
present. After 1:00 UT on Nov. 1, 9b and (pB evolve progressively 
away from their values found at the rear of the magnetic cloud, 
and their fluctuation level increases. This behavior is character- 
istic of the back region of a magnetic cloud that displays inter- 
mediate properties between a magnetic cloud and the solar wind 
dPasso et al.| f2006, 2007). 

According to Rich ardson & Cane| ( |2010| l the ICME ends 
when /3p shows a sharp increase. Nevertheless, even though /3p 
increases, it remains less than 0.4 and the ratio Tobs/Tg^p is still 
less than 0.5 almost until the end of Nov. 2 (Figure |2]i. In addi- 
tion, both 9b and (ps still display similar mean values as at the 
rear of the magnetic cloud. Thus, the back region extends until 
the end of Nov. 2 (when Tobs ~ T^exp and /3p ~ 0.4). We conclude 
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that the energetic particles of GLE 67 propagate up to the Earth 
in the extended back region of the previous magnetic cloud. 

Following a similar reasoning, we identified the magnetic 
structures of the IMF for the other GLEs (the analysis is re- 
ported in Section A. 2 of the Appendix). We conclude that rela- 



tivistic paiticles related to the GLE 61 propagate in the sheath of 
an ICME (Section A. 2. 2 1. The GLE 66 occurs inside a magnetic 
cloud (previously studied by Mandrini et al. |2007 1. It is worth to 
notice here that this magnetic cloud is related to the CME ejected 
from the Sun on Oct. 28, 2003 at 11:30 UT during the solar 
eruption producing the GLE 65. Therefore, we conjecture that 
the relativistic particles produced during the flare/CME event on 
29 Oct. 2003 have been injected in the footpoints of the CME 
launched on 28 Oct. 2003 (Section |A34|. For the GLE 59 and 
GLE 65, the temporal evolution of magnetic field and plasma 
parameters of the interplanetary medium displays intermediate 
properties between solar wind and magnetic cloud. A detailed 
analysis of all variables in a 4 days time interval around the GLE 
suggests that the particles related to GLE 59 and GLE 65 propa- 
gate in the back region (see Section A.2.1| and A. 2. 3 i. 

The four GLEs (59, 61, 65, 66) are all poorly connected to 
Earth by a nominal Parker spiral (Table [TJ. The solar eruptions 
associated to GLE 59, GLE 65 and GLE 66 are located near the 
central meridian, while the parent active region of GLE 61 is 
probably located behind the solar west limb. Thereby, the con- 
nection between the active region and Earth can not be ensured 
through the Parker spiral field line. According to our analysis, we 
concluded that particles propagate in transient magnetic struc- 
tures (ICME, magnetic cloud and back region). Since the coronal 
roots of such transient structures extend over a large longitudinal 
range, they can provide a magnetic path connecting the active re- 
gion and the Earth. However, this does not imply that all GLEs 
associated with transient magnetic structures should be poorly 
connected. Indeed, the broad longitudinal extent of an ICME in- 
cludes the well connected cases. An example is GLE 67 which is 
well connected by a nominal Parker spiral, while the high energy 
particles are traveling in a magnetic cloud back region. 



2.3.3. Disturbed solar wind 
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Fig. 3. Interplanetary data around GLE 60: example of a 
GLE present in a S W disturbed by a previous ICME. Same draw- 
ing convention as Figure [T] 



The identification of the interplanetary magnetic structure of the 
GLE 60 and GLE 69 indicates that particles propagate in an in- 
terplanetary medium formed by a disturbed solar wind. 

We present in detail below the case of the GLE 60, occur- 
ring on Apr. 15, 2001 (Figure [3]). The first relativistic particles 
impact the Earth at 13:55 UT. Between 10:00 and 16:00 UT, 
even though the magnitude of the magnetic field and the plasma 
parameters are close to typical solar wind values, the magnetic 
field orientation (Ob ~ -45° and <pB ~ 0°) does not correspond 
to the Parker spiral. However, the fluctuations of the magnetic 
field around fixed values of latitude and longitude are more char- 
acteristic of a magnetic stiTicture associated with a quiet solar 
wind than with an ICME. From the ICME list of Richardson 
|& Cane| ( |2010| l, the last strong magnetic perturbation was an 
ICME, starting the Apr 13 at 07:34 UT and ending on Apr 14 at 
12:00 UT, while the next ICME starts about 3 hours after the rel- 
ativistic particles reached Earth. Then, we suggest that particles 
propagate in a solar wind strongly disturbed by the preceding 
ICME. 

Similarly, the magnetic structure in which the relativistic 
particles of GLE 69 propagate is identified to be also a so- 
lar wind magnetic topology disturbed by a previous ICME 
(Section|A3]T). 



Judging from their longitudes, GLE 60 (W 85) and GLE 69 
(W 55) are, respectively, poorly and well connected to the Earth 
by the Parker spiral (Table[T]). However, the apparently good con- 
nection of the active region related to the GLE 69 does not imply 
that the magnetic field connecting the Sun to the Earth is a Parker 
spiral. The anisotropy of this GLE, with the strongest and fastest 
rise seen by neutron monitors in Antarctica, suggests a strong 
out-of-the ecliptic component of the magnetic field, which is ac- 



tually observed (Section A. 3.1 



2.4. Summary 

Using ACE magnetic field and plasma parameter measurements, 
we identified the structure of the interplanetary magnetic field. 
We show that, only in 3 events out of 10, relativistic particles 
possibly propagate in a solar wind structured by the Parker spi- 
ral, while in the 7 other events, relativistic particles propagate 
in a transient IMF (ICME, magnetic cloud, back region or dis- 
turbed solar wind). Among these 7 events, 5 have a solar source 
region significantly far (more than 30°) from the theoretical lo- 
cation of the footpoint of the Parker spiral at the solar source 
surface. 
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3. The interplanetary length 

Previous studies of relativistic particles, measured at Earth by 
neutron monitors, assume that energetic particles travel along the 
Earth-connected Parker spiral field line. This assumption implies 
that particles that are not scattered travel roughly 1.1-1.3 AU, de- 
pending on the solar wind speed. In fact, from our above results, 
most GLEs do not propagate in quiescent solar wind. Thus, we 
expect that the interplanetary length travelled by energetic parti- 
cles should difl'er from the one of the Parker spiral. 



3. 1 . Velocity dispersion analysis using the rising phase 
method (RPM) 

A commonly used estimate of interplanetary travel paths of en- 
ergetic particles is the velocity dispersion analysis (VDA) of the 
time when the first particles arrive at the detector The method is 
based on the assumption that the first energetic particles are si- 
multaneously released from a small acceleration site in the solar 
corona and propagate scatter-free in the interplanetary medium. 
Therefore, the less energetic particles arrive at the Earth later 
than the most energetic ones, implying a velocity-dependent 
time dispersion of the first particle detection. Plotting the arrival 
time of the first particles measured at Earth as a function of the 
inverse of particle speed (1/v), one predicts a linear relationship 
whose slope is a measure of the interplanetary distance travelled 
by the energetic particles, and whose intercept with the 1/v - 
axis is the solar release time. 

The major practical problem of the method is the definition 
of arrival times of particles at the detector The onset time is 
commonly defined as the instant when the intensity exceeds a 
given level above background. But this value is strongly affected 
by fluctuations in the individual energy channels, which may in- 
troduce a distortion and produce a large error both in the slope 
(the length) and the intercept (the release time) of the result- 
ing plot. This is what we experienced when applying different 
such methods to the events under study. We therefore propose a 
new method, called the rising phase method (RPM). The RPM 
is based on the same principle as the VDA, but it compares well- 
defined reference times during the rise phase of the time profiles 
at different energies, and thereby provides an estimate of the 
interplanetary length that is much less sensitive to background 
fluctuations than the classical VDA. Practically we proceed as 
follows: 

1. The intensity profile in each energy channel is divided by 
the background intensity averaged over one hour before any 
evident rise of the signal. Since the signal has a broad range 
of variation we take the logarithm of this ratio, making the 
background level zero. 

2. Depending on the energy spectrum of the event, the different 
energy channels reach different maximum values. In order 
to compensate this hardness effect, we normalize the above 
logarithm of fractional intensity by its maximum found just 
after the rise phase. This normalization is realized in a time 
interval typically between 10 and 20 min after the rise phase, 
with a time-shift function of the channel energy in order to 
compensate approximately the differential arrival time in the 
energy channels. Thus, the logarithm of this normalized in- 
tensity increases from to 1 in the rising phase of the signal. 
Three examples are shown in the top row of Figure]?] 

3. The rising phase between the levels 0.2 and 0.8 is fitted by 
a straight line (assuming an exponential rise of intensity). 
The fit reduces the effects of the fluctuations superposed on 



the general increase of intensity, as seen in several curves 
in Figure ]4] panels a & c. We refer to the instant when the 
fitted straight line is at level 0.5 as the rise time at the cor- 
responding energy. This corresponds to the maximum cross- 
correlation of the rising profiles. 
4. Finally, we plot the rising parts of intensity profiles as a func- 
tion of observing time minus the rise time, and control vi- 
sually how well the profiles in the different energy ranges 
superpose. 

As the classical VDA, the RPM provides a set of times frise 
which depend on the energy of the channels, hence on the veloc- 
ity of the detected energetic particles. Then, the slope of a linear 
fit of fiise versus 1/v provides an estimation of the interplanetary 
length D (bottom row of Figure]4|i. The estimation of D depends 
only on the relative timing of the channels. 

The linear fit also provides an estimation of the solar release 
time, fsRT of the particles (supposed to be independent of energy) 
as the intercept of the straight line with the 1/v = axis. We have 
to keep in mind that this is an upper limit, because it refers to a 
time during the rise of the intensity profile. We estimate the error 
on D and fsRx with the standard error of the fitted parameters. 

We believe that the use of the rise times gives a more reli- 
able determination of the travel path than using the onset times 
mainly for three reasons. As akeady stated, the method does not 
depend on the fluctuating background, since it uses the part of 
the signal well above it. Moreover, the fluctuations of the back- 
ground and of the signal should have a minor influence since 
the method cumulates the information over a finite time inter- 
val (most of the rising part of the signal). Finally, the rise time 
is computed in the steepest part of the rising phase of the time 
profile, where the timing is the most accurate. These three prop- 
erties contrast with the onset time method which relies on a local 
measure above the fluctuating background where the signal just 
starts to increase. 

However, the time dispersion determined by the RPM may 
be affected by systematic biases, such as scattering during 
the propagation of particles from the Sun to Earth. Since the 



mean free path increases with increasing particle energy ( Droge 



|2000 1, the lower energy particles should arrive with an extra de- 
lay compared to the arrival time that one expects. This extra de- 
lay at low energy may lead to an increase of the length obtained 
by the RPM. The dependance of the RPM's results in respect to 
the particle scattering is hard to assess. However, we expect that 
one of the following three situations should apply: (1) scattering 
of the first arriving particles is negligible, and the slope of the 
onset time vs 1/v plot gives the geometrical path, e.g. GLE 63 & 
GLE 64; (2) scattering dominates the propagation to the extent 
that the slope is mostly determined by the energy dependence of 
the scattering mean free path, e.g. GLE 70 ; (3) both the geomet- 
rical path and scattering shape the plot. Thereby, the RPM gives 
at least the apparent length travelled energetic particles, but this 
apparent length may be a mix of the geometrical length and the 
additional length created by the scattering. 

3.2. Application to the selected events 
3.2.1 . Data set and event selection 

We perform the velocity dispersion analysis, based on the RPM, 
using measurements of the High Energy Detector (HED) of the 
Energetic and Relativistic Nuclei Experiment (ERNE) instru- 
ment ( jTorsti et aTj |1995| ) on board SoHO, located at the LI 
Lagrangian point. The energy range covered by ERNE/HED 
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Fig. 4. Example of the RPM and VDA analysis, top row: normalized time profile of the logarithm of the intensity during the lising 



phase without any time shift, for GLEs 64, 67 and 70, respectively on panels a, b and c (see Section 3.1 1. The energy range (in MeV) 



of channels is added on the panel a. Bottom row, panels d, e and f: plots of rise and onset times in function of the inverse of the 



proton velocity. The dark grey (or blue) points correspond to the rise times determined by RPM apply to ERNE data (Section 3.1 



and the related linear fit is displayed by a dashed line (dark grey or blue). The linear fit related to the data set including the NM 
point (labeled and plotted in black or red), in addition to the previous ERNE points, is represented by a continuous line (black or 
red). The medium grey (or green) points correspond to the onset times of ERNE data and the related linear fit is displayed with a 
dashed line and the same color The linear fit with the NM point added is shown with a continuous line and light grey (or pink). A 
color version is available in the electronic version. 



extends from 13 MeV to 130 MeV. Data are provided in 10 
customized energy channels. The 5 highest energy channels, 
between 40 - 130 MeV, are usually distorted during GLEs 
(R. Valtonen, private communication): in order to privilege the 
detection of heavy ions, the onboard software raises the thresh- 
old for the detection of protons and He ions during large events. 
As a result, the proton intensities between 40 - 130 MeV may 
display recurrent indentations. We correct this artifact in order to 
recover a regular time profile to which the RPM can be applied. 
Identifying the start and end times of the indentations by the 
sudden change of intensity, we apply a correcting factor to each 
interval by using the measured values obtained just before and 
after the jump. These coiTections have a negligible effect on the 
fitted parameters. Still, they allow to use higher energy channels 
which are closer to the NM data. 

Since it is not clear a priori if the relativistic protons mea- 
sured by the neutron monitors are just the high-energy exten- 
sion of the SEP spectrum measured by ERNE or a distinct pop- 
ulation, we perform separate velocity dispersion analyses with 
and without the data point of the GLE as determined by Moraal 
& McCracken| ( |20TTl l. We use the mean energy of each ERNE 
channel to compute the velocity of energetic protons. The RPM 
can not be applied to neutron monitor measurements, which re- 
spond to the integral proton spectrum above the local rigidity 
cutoff, as determined by the geomagnetic field configuration. 

Only seven of the ten GLEs have suitable ERNE observa- 
tions for the RPM: there is no data for the GLE 59 and for the 
GLE 66, energetic particles of the SEP are smothered in a larger 
energetic particles flux probably accelerated at the bow shock of 
the magnetic cloud. Finally, we can neither use the ERNE data 



Table 2. Path lengths and release times using RPM 



1 2 3 

# GLE Interplanetary length 



4 5 
Solar release time 





- NM, rise 


-1- NM, rise 


- NM 


rise 


+ NM, rise 


60 


1.50±0.16 


1.55 ± 


0.10 


13:45 


±05 


13:43±03 


61 


1.61 ±0.12 


1.63 ± 


0.09 


02:23 


±04 


02:23±03 


63 


1.27 ±0.06 


1.39 ± 


0.07 


05:35 


±02 


05:3 1±02 


64 


1.31 ±0.09 


1.32 ± 


0.06 


01:12 


±03 


01:12±02 


65 


1.59 ±0.16 


1.89 ± 


0.17 


11:15 


±05 


11:05±05 


67 


2.67 ±0.15 


2.55 ± 


0.09 


17:02 


±05 


17:07±03 


70 


2.87 ±0.13 


2.56 ± 


0.16 


02:11 


±04 


02:20±05 



Note : Effective interplanetary length (AU) travelled by protons and the 
associated solar release time (UT). The error bars correspond to the 
standard deviation of the linear fit. Col. 1: the label of the GLEs. 
Col. 2: length computed from the rise times and the mean energy 
of ERNE channels - without the NM's data. Col. 3: as Col. 2 but 
with the NM's point. Col. 4: the solar release time of energetic pro- 
tons from the RPM applied to only to ERNE data (as Col. 2, see 
Figure|4]l. Col. J: as Col. 4 but with the NM's point. 



of the GLE 69 because the signal decreases and no bump has 
been detected in the time interval during which we expect the 
arrival of energetic particles related to the GLE. 



3.2.2. Results 

Figure|4]displays the results of the RPM for three GLEs selected 
to illustrate the three kinds of interplanetary magnetic structures 
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Table 3. Path lengths and release times using the classical VDA 



1 

#GLE 



3 

Interplanetary length 





- NM, onset 


+ NM, onset 


Reames 


60 


0.93 ± 0.26 


0.85 ±0.17 


1.59 ±0.01 


61 


1.60 ±0.49 


1.34 ±0.38 


1.80 ±0.10 


63 


0.78 ± 0.07 


0.76 ± 0.05 


1.64 ±0.06 


64 


0.83 ±0.17 


0.82 ±0.13 


2.16 ±0.05 


65 


3.12 ±0.44 


2.41 ± 0.43 


1.38 ±0.03 


67 


2.65 ± 0.23 


2.27 ±0.18 


2.01 ±0.04 


70 


2.47 ± 0.21 


1.91 ±0.29 


2.81 ±0.02 



Note : Effective interplanetary length (AU) travelled by protons and 
the associated solar release time (UT). Col. 1: the label of the GLEs. 
Col. 2: length computed from the onset times and the maximal en- 
ergy of ERNE channels - without the NM's point. Col. 3: as Col. 2 
but with the NM's point. The error bars of Col. 2 and 3 correspond 
to the standard deviation of the linear fit, increased by the error 
on the determination of the onset time (due to background fluctua- 
tions). Col. 4: length from |Reamesl ( |2009[ >. 



previously identified. Even though the intensity profiles display 
significant fluctuations (top panels of Figure |4]), the RPM suc- 
ceeds to derive the rise time of energetic particle fluxes in each 
energy channel. When these rise times are plotted as a function 
of 1 / V (bottom panels of Figure]?]), they display a linear relation. 

We repoit in Table]2]the interplanetary lengths resulting from 
our RPM with and without the NM data, respectively. We found 
that for 5 out of 7 events, the difference between the two lengths 
computed with and without the neutron monitor data point is not 
larger than ^ 0.05 AU, so well inside the error bars. 

However, we emphasize that our RPM considers an instant 
near half maximum during the lise of the ERNE time profiles 
as the onset time, whereas the neutron monitor onset time is 
the instant when the first neutron monitor signal was detected. 
These times are not directly comparable. A number of the events 
of Table |2] for which ERNE measured well-defined onset time 
profiles, are associated to weak GLEs, where the time when the 
signal starts to exceed background is probably a rough upper 
estimate of the actual start, and actually also designates a time 
during the rise of the GLE. It is likely that the systematic error of 
the neutron monitor onset time compensates to some extent the 
different timing definitions of the two data sets. Therefore we 
note that the timing is consistent, but we can not use this finding 
to draw firm physical conclusions on the simultaneity of proton 
acceleration from tens of MeV to GeV energies. 

A discrepancy between the results with and without NM 
data is found in GLEs 65 and 70. They display a difference 
comparable to the error bar sizes (Table |2] and Figure ]5]). For 
GLE 70 the results with and without NM has a length differ- 
ence of ^ 0.31 AU. Irrespective of the exact numerical value, 
both analyses lead to an interplanetary path that is significantly 
longer than the Parker spkal for GLE 70. 

In the case of the GLE 65, we adapt the RPM. Indeed, the 
normalized logarithm of intensity has a temporal evolution sig- 
nificantly dependent of the energy channel and so has not a sim- 
ple temporal shift as for other GLEs. We tested the usual nor- 
malization method, but it appears that the rising phases of the 
time profiles in different energy channels are not parallel. This 
difference in the time evolution of the channels indicates an evo- 
lution of the hardness of the spectrum during the rise phase. 
Indeed, one notices that a pre-increase appeared in the ERNE 
time profile in the high energy channels. According to |Trottet| 



|et al.|p008] l, the acceleration of particles during this solar erup- 
tion displays a complex temporal structure, suggesting that sev- 
eral episodes of acceleration occur. Thereby, the pre-increase can 
result from a distinct acceleration episode than the one acceler- 
ating the particles of the main particle flux. 

In order to remove this effect for GLE 65, we use an earlier 
normalization, during the rising phase, at a time chosen such as 
to obtain rise time profiles in different energy channels that are 
paraUel. We also find that the rise times are more aligned with 
the neutron monitor data and that the lengths computed by the 
RPM with and without the NM point differ only by ^ 0.3 AU. 

3.2.3. Comparison with results of classical VDA 

We compare in this section the results of the RPM with the VDA 
method based on the estimation of the onset time, fonset, of the 
first particles in each channel. Since the arrival time of the first 
protons is masked by the fluctuations of the background, we ex- 



trapolate the linear fit of the rise phase (see Section 3. 1 1 up to the 
mean level of the background. This defines the onset time, fonset- 
The uncertainty on fonset is partly due to the error on the linear 
fit of the rising phase (as for fiise), but it is also increased by 
the extrapolation toward the background level and by the back- 
ground's fluctuations. Considering the standard deviation of the 
background, this extra uncertainty on fonset is typically of 3 min. 

Since the first particles are expected to have the highest en- 
ergy, we compute 1/v from the maximum energy of each chan- 
nel. Using the maximum energy, fimax, compared to the mean 
one, /imean introduccs a systematic difference in the estimated D 
of a factor 1 -H A£'/(2£'mean) ~ L06 to the first order in A/i/ZSmean 

(where ^E = E^^y, - E,^i„ and Emean = (^max + £'min)/2). 

We find that the results with fonset are typically much more 
fluctuating and incoherent with the neutron monitor data than 
the results with fiise (Figure]4|. Moreover, the length found with 
the onset times is unrealistically short (below 1 AU) for some 
GLEs (Table|3]l. 

In Table ] 3] w e also report the interplanetary length com- 
puted by Rea mes| ( |2009| l with a velocity dispersion analysis of 
the onset time of protons and heavy ions. These results have 
some large differences with our VDA results (Table]3]l. However, 
Reames results are broadly consistent with our rise time results 
(Table |2|, except mainly for three events: GLEs 63, 64 and 67. 
Even though the VDA has been applied to the heavy ion data 
(from WIND/EPACT/LEMT) in order to get a lower background 
level, the fluctuations of the signal before and during the rising 
phase, see e.g. Figure 3 in Reames (2009 1, and the temporal res- 
olution (5-10 min) introduce limits to the accuracy of the esti- 
mated lengths. 

In summary, we conclude that RPM gives a more meaning- 
ful evaluation of the interplanetary length travelled by energetic 
particles than the VDA, whatever the data set is. This conclusion 
is further strengthened in the next sub-section. 



3.3. Comparison of the RPM results with the identification of 
the IMF structure 

3.3.1. IMF structure and interplanetary path length 

From the comparison of the interplanetary path length inferred 
from the RPM (Table ]2]| and the interplanetary structure iden- 
tified in Table [T] a clear relation can be established between 
the magnetic field along which energetic protons propagate and 
the distance they travelled. Figure ]5] summarizes the interplan- 
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Fig. 5. Summaiy of the estimated length by the RPM, with (cross 
symbol) and without NM (diamond symbol), in relation with the 
magnetic topology. 



etary lengths computed without and with NMs (Col. 2 and 3 in 
Table [2| as a function of the IMF topology (Col. 8 in Table [T}. 

The protons in the rising phase of GLE 63 and GLE 64 that 
are associated to a Parker spiral magnetic structure, travelled ap- 
proximately 1.3 ± 0.1 AU. These interplanetary lengths are con- 
sistent with the Parker spiral magnetic topology that has a theo- 
retical length of * 1.15 AU for a slow solar wind. These results 
differ from Reames ( 2009 | l who found a length of » 1 .6 AU and 
2.2 AU (Table [3]) for GLE 63 and 64, respectively. Within the 
error bars (+0.06), these lengths are compatible with a Parker 
spiral only if an extra delay is present, so if a physical mecha- 
nism can provide it independently of the proton energy. 

By contrast of GLE 63 and 64, the interplanetary lengths 
travelled by protons producing GLEs 60, 61, 65, and 67, are 
longer than 1 .5 AU (see Table|2]l. Such lengths are fully compat- 
ible with the magnetic topology of an IMF created by transient 
magnetic structures, such as an ICME, magnetic cloud, back re- 
gion and even a disturbed solar wind. 

The two lengths obtained with and without NMs, for GLE 65 
(see Table |2|, are both consistent with the back region of a 
magnetic cloud (because a back region is formed by reconnec- 
tion between the solar wind and the flux rope, see the end of 
Section [2T|). However, the longer length is supported by the re- 
sults from [Miroshnichenko et al. ( 2005| l. Studying this specific 
GLE, they showed that the delay of the proton arrival should be 
due to the perturbation of the IMF and that protons should travel 
roughly 2 AU. In addition, they estimated that this length is con- 
sistent with a simultaneous release of relativistic protons with 
relativistic neutrons, estimated at 10:56 UT from the NM data 
(tsRT - Inm - D/v). From our RPM analysis, and for linear re- 
gression with and without NMs, the proton release occurs later 
(Table |2]l. However, the solar release time may have been over- 
estimated because the pre-increase detected in the high energy 
channels of ERNE has been ignored in the determination of the 
time dispersion (Section [3.2.2| i. 
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Fig. 6. Comparison of the solar release time of relativistic pro- 
tons obtained from the RPM with (cross symbol) and without 
NM (diamond symbol), with time interval of electrons injection 
as deduced from type III (grey shaded bars). The onset time of 
type III bursts is set as the origin of the vertical axis. 



The last relativistic event of the solar cycle 23 (GLE 70) 
displays also a discrepancy between the magnetic structure and 
the interplanetary length. Indeed, the length travelled by protons 
is a; 2.6-2.9 AU, while the IMF displays characteristics of the 
Parker spiral. We note that this is the only GLE of our sam- 
ple that is detected while the Earth is in the fast solar wind. It 
is well known from Ulysses observations in the fast solar wind 
that SEP are more strongly scattered there than in the slow wind 
(Sanderson 2004). Since the mean free path of the particles in- 
creases with increasing energy, particles of lower energy will 
be more strongly delayed than if they propagated scatter-free. 
This may contribute to an overestimation of the interplanetary 
path length. Moreover, the solar source of the fast solar wind is 
typically further from active regions than the edge of the open 
flux region, so from the source region of a GLE. Then, when 
energetic particles are detected in the fast solar wind they are 
expected to have been transported a long way from their acceler- 
ation site. This extra delay could contribute to the long apparent 
length found for GLE 70. Such effect is expected to be smaller 
in ICME and slow solar wind cases as the guiding magnetic field 
is located closer to the GLE solar source. 

Except for this complex GLE 70, the interplanetary length is 
consistent with the magnetic structure observed during the parti- 
cle event. We conclude that the propagation of energetic protons 
occurs in an IMF specific to each GLE, and frequently different 
from the Parker spiral often used. Moreover, this shows that the 
geometry of the magnetic structure affects the length travelled 
by particles from the Sun to the Earth. 



3.3.2. When are protons injected into tine IMF? 

The propagation of non-thermal electrons in the interplanetary 
medium generates radio emission in the decametric wavelength 
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range with a characteristic drift towards lower frequencies as the 
electrons propagate away from the Sun, towards lower ambient 
densities, along open magnetic field lines. These are type 111 
bursts. Large groups of type III bursts accompany most major 
SEP events ( Cane et al\ |2002| l. At 14 MHz, it is generally as- 
sumed that the electromagnetic waves are emitted at ~ 2 - 3 R© 
from the solar surface. The timing of the type III bursts at this 
frequency gives an indication of the time of electron injection 
into the high corona. Although the electrons are accelerated 
lower in the corona, the travel time from the acceleration site 
to the 14 MHz source is only a few tens of seconds, so this cor- 
rection of time is negligible. 

We computed the time interval of electron injection with the 
WIND/ WAVES [Bougeret et al.| ( [T995) data at 14 MHz and we 
compared them to the solar release time. The Figure |6] synthe- 
sizes the temporal relation of the injection of the non-thermal 
electrons and the relativistic protons for the 7 studied events 
(Section [TZ2| . 

For 6 of the 7 events studied, the solar release time of ener- 
getic protons is enclosed in the time interval of electron injection 
(Figure[6]l. For GLE 67 the solar release of energetic protons pre- 
cedes that of electrons by 3-5 min. This is within the uncertainty 
of the method (see Table|2]l. We therefore conclude that the first 
energetic protons during these events are released during the in- 
terval of the electron injection traced by the type III bursts. 

This timing consistency is an additional argument support- 
ing our results on the magnetic topology of the interplanetary 
medium and its role for the propagation of energetic particles. 
Essentially, it gives a second proxy confirming that the inter- 
planetary length has been well evaluated. 



4. Conclusion and Discussion 

We performed independent studies based on different methods, 
in order to define and constrain the characteristics of the inter- 
planetary magnetic field during GLEs. The main results are sum- 
marized below. 

First, in most of the events, relativistic protons propagate in 
transient magnetic structures and not along the Parker spiral field 
lines. 

Second, the magnetic structure of the IMF is consistent with 
the location of the parent active region of the GLE: when the 
in situ measurements show a Parker spiral, the parent active re- 
gion is within ±30° of the nominal footpoint, whereas it may be 
farther away when the IMF has a transient configuration. 

Third, the effective interplanetary length computed through 
our velocity dispersion analysis (rising phase method, RPM) for 
energetic and relativistic protons is in agreement with the mag- 
netic topology of the IMF deduced from in situ measurements. 
This length travelled is 1.3 AU when the protons travel in the 
slow solar wind while it is in the range 1.5-2.6 AU for transient 
magnetic structures. For the only GLE having protons traveling 
in the fast solar wind, the effective length is much longer than ex- 
pected (2.6-2.9 AU rather than 1.1 AU). This may indicate that 
an extra mechanism delayed the protons the more strongly, the 
lower their energy. Pitch angle scattering is a possible candidate, 
especially since it is expected to be stronger in the fast than in 
the slow solar wind. 

Fourth, our results are consistent with the idea that the first 
protons producing SEPs and GLEs are accelerated/released si- 
multaneously into the interplanetary medium, although they do 
not demonstrate this. The weakness of neutron monitor signa- 
tures in a number of events where the RPM could be applied to 



ERNE data precludes firm conclusions. The analysis also sup- 
ports the idea that the first arriving energetic protons are not sig- 
nificantly scattered during their travel from the Sun to the Earth. 

Fifth, according to the solar release time, energetic protons 
are injected during the injection of non-thermal electron beams 
producing type III bursts, supporting the above results on the im- 
pact of the interplanetary magnetic structure on the proton trans- 
port. 

The common release of relativistic protons and type Ill- 
related electrons has already been suggested in a detailed study 
of GLE 69 ( |Masson et al.||2009] l. The event was not studied here, 
because the interplanetary magnetic field measurements were 
unreliable. But the result is fully consistent with the conclusions 
drawn from the seven GLEs of the present sample. Our observa- 
tional study suggests that the common release of type III bursts 
emitting electrons and the first escaping energetic and relativistic 
protons is a general property of GLEs. 

In addition, energetic particles related to active region poorly 
connected to the Earth through the Parker spiral are may con- 
nected through a transient magnetic structure. This result implies 
to reconsider the driven-shock acceleration being invoked to ex- 
plain the detection of particles at Eaith in poorly connected SEP 
events. 

Finally, energetic protons propagating in a transient mag- 
netic structure travel a longer distance than they should do in 
the Parker spiral case. This delay, due to the geometrical length 
of the interplanetary magnetic field, contributes to understand 
the quasi-systematic delay measured between the GLE and the 
first radiative signature of accelerated particles at the Sun (jChver 



et al. 1982 1. Moreover, this longer length travelled by relativis- 



tic protons modifies the relation between the coronal radiative 
signatures and the detection at Earth of these protons. As an 
example, Li et al.| (2009') assumed that relativistic protons of 
GLE 70 travelled 1.1 AU between the acceleration site and the 
Earth, leading them to conclude that particles have been flare- 
accelerated. However, our study clearly demonstrates that for 
GLE 70, one can not consider such a short path length of the 
interplanetary magnetic field. 
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Fig.A.l. Interplanetary data around GLE 64: it is located in a 
Parker like solar wind. Same drawing convention as Figure[T] 



Appendix A: Interplanetary magnetic structures for 
other GLEs 

This appendix provides the analysis of the interplanetary struc- 
tures associated with those GLEs which have not been detailed 
in the main text. 

A.1. Solar Wind 

A.1 .1 . GLE 64 on Aug. 24, 2002 

Relativistic protons of the GLE 64 arrive at Earth at 01 :24 UT on 
Aug. 24, 2002, during the decrease of B from 12 nT to 5 nT in 
This decrease could suggest that protons 



two days (Figure A. 1 



travelled in the wake of the previous ICME ending on Aug. 21, at 
14:00 UT ( [Richardson & Cane| [2010) 1. However, the orientation 
of the magnetic field 0b ~ 0°, cpB - -45°, the criteria on the 



temperature, T, 



exp 



Tobs, and /3p - 0.4, are consistent with the 



Parker spiral magnetic field. 



A.1.2. GLE 70 on Dec. 13, 2006 



The magnetic field does not show any increase of magnitude or 
magnetic rotation, between 00:00 UT on Dec. 12 and 24:00 UT 
on Dec. 13 (Figure |A.2|i. Indeed, the magnetic field is non- 
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00:00 08;00 16:00 00:00 08:00 16:00 00:00 
Start Time : 06/12/lS, 0000:00 UT 

Fig.A.2. Interplanetary data around GLE 70: it is located in a 
Parker like solar wind. Same drawing convention as Figure[T] 
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Start Time : 00/07/13,0000:00 UT 

Fig. A.3. Interplanetary data around GLE 59: it is located in the 
back region of a hot flux rope. Same drawing convention as 
Figure [T] 



coherent and fluctuates around 6b = ° and <pB = -45 °, 
with stronger fluctuations before the GLE and weaker after that. 
Moreover, the expected temperature is roughly equal to the ob- 
served temperature, and 0.4 ^ /3p ^ I. The magnetic and plasma 
properties strongly suggest that the Parker spiral is the mag- 
netic structure of the IMF, guiding the relativistic protons of the 
GLE 70, detected at the Earth at 02:50 UT on Dec. 13, 2006. 

A.2. ICME, magnetic cloud, or back region 
A.2.1. GLE 59 on Jul. 14, 2000 

The magnetic field on Jul. 14, 2000 between 02:00 UT and 
15:00 UT, enclosing the arrival time of the GLE 59 detected at 
10:31 UT, does not display any common features of the well 
defined interplanetary magnetic field (Parker spiral, ICME, or 



2000 at 13:00 UT up to Jul. 14 at 02:00 UT, which is typical of a 
magnetic cloud with its sheath in front. Meanwhile, Tobs > Texp 
indicates that this is a hot flux rope. Even being hot, it was previ- 



magnetic cloud. Figure A.3i. The magnetic field is weak {B ^ 
3-5 nT), while the magnetic vector orientation is highly fluc- 
tuating and so does not have the typical Og and cpg values of the 
Parker spiral. The plasma parameter evolution before the GLE 
also displays inconsistency: Tobs - 0.5rexp, like in an ICME, 
whereas 0.4 < /3p < I, more typical of the solar wind. 

Before this time period, the magnetic field magnitude in- 
creases and its magnetic vector rotates coherently from Jul. 13, 



ously classified as a magnetic cloud ( Huttunen-Heikinmaa et al. 
2005 ; Ri chardson & Cane} 2010[ l. These previous studies set the 
ICME end on Jul. 14, 2000 at 15:00 UT. The coherent rotation 
of the field, so the flux rope, ends at a!02:00 UT. Later on the 
magnetic field orientation has large fluctuations. Still, globally 
(pB progressively evolves from its value at the rear of the flux 
rope ((pB ~ 250 - 280°) to an outward sector value. Thus, we 
propose that the region between 02:00 UT and 15:00 UT is the 
magnetic cloud back region ((pB ~ 135°). The magnetic field and 
plasma properties are intermediate between the solar wind and 
magnetic cloud. Therefore, our analysis suggests that relativistic 
protons of GLE 59 arrived at Earth in the back region of a hot 
flux rope (magnetic cloud-like). 



A.2.2. GLE 61 on Apr. 18, 2001 

Relativistic protons on Apr. 18, 2001 (GLE 61), reach Earth at 
02:36 UT, during a strong increase of the magnetic field strength 
(B ^ 20 - 25 nT, Figure |A.4[). After this period, we identify an 
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Fig. A.4. Interplanetary data around GLE 61: it is located in the 
sheath of an ICME. Same drawing convention as Figure [T] 
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Fig.A.5. Interplanetary data around GLE 65: example of a 
GLE present in the back of a magnetic cloud. Same drawing 
convention as Figure [T] 



ICME (B ~ 10 nT, decreasing progressively to ~ 5 nT) from 
the Apr 18 at 12:00 UT to the Apr 20 at 11:00 UT. This tem- 
poral evolution of B is the signature of an ICME, preceded by 
its sheath. Our identification is in agreement with the boundaries 
determined by composition analysis ( Richardson & Cane,,2010 ). 
However, this ICME is a particular case where the proton tem- 
perature is higher than the expected one and therefore /3p ^ 1 
in the front part of the ICME (Figure A.4i. We conclude that 
the relativistic protons of GLE 61 traveled in the sheath of this 
ICME. 



A.2.3. GLE 65 on Oct. 28, 2003 

Figure [A3] displays the evolution of the interplanetary magnetic 
field and plasma parameters for the GLE 65. The increase and 
the temporal evolution of the magnitude of the magnetic field 
starting at 02:00 UT on Oct. 28, 2003, the high coherence of the 
magnetic field vector, the ratio of the expected and the observed 
temperatures approaching 2 and the l ow value of /3p <QA clearly 
indicate that ACE crosses an ICME. Richardso n & Cane| ( |2010| l 
determined its end at 09:00 UT on Oct. 28. They do not classify it 
as a magnetic cloud probably because the magnetic field rotation 
is only » 60°. This limited rotation indicates that the flux rope 



was crossed by the spacecraft with a large impact parameter (i.e. 
far from its axis). Then, we conclude that the in-situ data are 
compatible with a magnetic cloud. 

Around 09:00 UT the magnetic vector has a strong discon- 
tinuity, followed by many others until 16:00 UT on Oct. 28, as 
well as a progressive evolution of (pg- During this time interval, 
the observed temperature becomes about twice higher than the 
expected temperature. The evolution of the magnetic field and 
of plasma parameters after 09:00 UT on Oct. 28 are compati- 
ble with the mixed properties of a hot back region in the wake 
of the magnetic cloud (Dasso et al. 2006 20071). However, be- 
cause the impact parameter of the magnetic cloud is large, the 
back region is also not observed in favorable conditions for a 
clear identification. Finally, relativistic protons are detected at 
Earth at 11:12 UT, and we thereby conclude that they propagate 
in this hot back region. 

A.2.4. GLE 66 on Oct. 29, 2003 

The interplanetary magnetic structure on Oct. 29, 2003 has al- 



ready been identified as a magnetic cloud (Mandrini et al. 



2007| ). The magnetic cloud starts at 1 1 :25 UT on 29 Oct. anc 
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:00 08:00 12:00 16:00 20:00 00:00 04:00 
Start Time : 03/10/29,0400:00 UT 

Fig. A.6. Interplanetary data around GLE 66: it is located inside 
a magnetic cloud. Same drawing convention as Figure[T] Plasma 
data are not available for this event. 



ends at 02:00 UT on 30 Oct. The relativistic particles of the 
GLE reached Earth at 21:01 UT on Oct. 29, then they propagate 



in this well-identified magnetic cloud (Figure A. 6 1. In addition, 
this magnetic cloud is related to the CME ejected from the Sun 
the 28 Oct. 2003 at 1 1 :30 UT during the solar eruption producing 
the GLE 65. Therefore, the relativistic protons produced during 
the flare/CME event the 29 Oct. 2003 have been injected in the 
foot-points of the CME launched the 28 Oct. 2003. 

A.3. Disturbed solar wind 
A.3.1 . GLE 69 on Jan. 20, 2005 

Before the GLE 69 on Jan. 20, 2005 at 6:49 UT, the interplane- 
tary magnetic structure is identified as an ICME, starting on Jan. 
18, at 23:00 U T and ending on Jan. 20, at 03:00 UT ( [Richardson 
& Cane 2010 1. During this time interval, the magnetic field has a 
significant strength (6-18 nT) with small fluctuations, and with a 
weak coherent rotation. Moreover, Tobs < 0.5rexp and /3p - 0.02. 
So all the characteristics defining a magnetic cloud, except a 
large rotation of the field, are present. This is a magnetic cloud- 
like case (Figure A. 7 1. 

After 03:00 UT, the magnetic field and plasma parameters 
behavior do not correspond to the Parker spiral case as the 
magnetic field is nearly orthogonal to the ecliptic plane and (ps 
is variable (so not characteristic of a solar wind sector). The 
large rotation of the magnetic field around 03:00 UT would 
be an extreme case for the start of a back region, difficult to 
explain by the reconnection of magnetic cloud with a solar wind 
magnetic field, then we cannot identify it to a back region. Still, 
the field nearly orthogonal to the ecliptic plane, the increase of 
the temperature and the Pp are consistent with a disturbed solar 
wind located in the wake of the preceding ICME. Even though 
any clear magnetic structure can be identify, one can at least 
certainly conclude that relativistic protons related to the GLE 69 
propagate in a disturbed solar wind. 
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Fig.A.7. Interplanetary data around GLE 69: it is located in a 
disturbed solar wind. Same drawing convention as Figure[T] 



Next page: Color version of one figure (for the elec- 
tronic version). 
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Fig. 4. 



Example of the RPM and VDA analysis, top row. normalized time profile of the logarithm of the intensity during the lising 

The energy range (in MeV) 



phase without any time shift, for GLEs 64, 67 and 70, respectively on panels a, b and c (see Section 3.1 



of channels is added on the panel a. Bottom row, panels d, e and f. plots of rise and onset times in function of the inverse of the 



proton velocity. The dark grey (or blue) points correspond to the rise times determined by RPM apply to ERNE data (Section 3.1 
and the related linear fit is displayed by a dashed line (dark grey or blue). The linear fit related to the data set including the NM 
point (labeled and plotted in black or red), in addition to the previous ERNE points, is represented by a continuous line (black or 
red). The medium grey (or green) points correspond to the onset times of ERNE data and the related linear fit is displayed with a 
dashed line and the same color The linear fit with the NM point added is shown with a continuous line and light grey (or pink). A 
color version is available in the electronic version. 



